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ABSTRACT 

In order to understand the physical mechanisms at work during the formation of 
massive early- type galaxies, we performed six zoomed hydrodynamical cosmological 
simulations of halos in the mass range 4.3 x 10 12 ^ M v i r ^ 8.0 x 10 13 M at z = 0, 
using the Adaptive Mesh Refinement code RAMSES. These simulations explore the 
role of Active Galactic Nuclei (AGN), through jets powered by the accretion onto su- 
per massive black holes on the formation of massive elliptical galaxies. In the absence 
of AGN feedback, large amounts of stars accumulate in the central galaxies to form 
overly massive, blue, compact and rotation-dominated galaxies. Powerful AGN jets 
transform the central galaxies into red extended and dispersion-dominated galaxies. 
This morphological transformation of disc galaxies into elliptical galaxies is driven by 
the efficient quenching of the in situ star formation due to AGN feedback, which trans- 
form these galaxies into systems built up by accretion. For galaxies mainly formed by 
accretion, the proportion of stars deposited farther away from the centre increases, 
and galaxies have larger sizes. The accretion is also directly responsible for randomis- 
ing the stellar orbits, increasing the amount of dispersion over rotation of stars as a 
function of time. Finally, we find that our galaxies simulated with AGN feedback bet- 
ter match the observed scaling laws, such as the size-mass, velocity dispersion-mass, 
fundamental plane relations, and slope of the total density profiles at z ~ 0, from 
dynamical and strong lensing constraints. 

Key words: galaxies: formation - galaxies: elliptical and lenticular, cD - galaxies: 
kinematics and dynamics - galaxies: active - galaxies: jets - methods: numerical 



1 INTRODUCTION 

Early-type galaxies (ETG) are amongst the most massive 
galaxies observed in our Universe. They are red, almost 
dead, and ellipsoidal-shaped objects lying in the densest re- 
gions of the Universe at redshift zero, in groups and clusters 
of galaxies. 

In the standard ACDM paradigm, progenitors of mas- 
sive collapsed structures such as groups and clusters of 
galaxies form at high redshift by the accretion of cold 
filamentary gas. At low redshift, the gas cooling rate at 
the virial radius becomes insufficient to efficiently evac- 
uate the energy of infall transformed into heat due to 
shocks, and leads to the formation of large, diffuse and 
hot quasi- spherical regions of gas around bright galax- 
ies (Birnboim & Dekel 2003; Keres et al. 2005; Ocvirk et al. 
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2008; Dekel et al. 2009). However, in the absence of any 
strong source of feedback from galaxies, theory fails to 
produce this population of red massive galaxies (e.g. 
Bower et al. 2006), as the gas cooling flows proceed unin- 
terrupted in the cores of halos. This leads to overly massive 
and actively star- forming galaxies at redshift zero, in stark 
contradiction with observations. 

Supermassive black holes (BHs) are commonly observed 
in the centres of galaxies with bulges. The BH mass Mbh 
scales with its host bulge stellar mass Mb, and stellar ve- 
locity dispersion ah (Magorrian et al. 1998; Tremaine et al. 
2002; Haring & Rix 2004). This has led several authors to 
suggest that a self-regulation process is taking place between 
BHs and their host galaxy (Silk & Rees 1998; King 2003; 
Wyithe & Loeb 2003), as accretion onto supermassive BHs 
can release tremendous amounts of energy that can poten- 
tially drive large-scale outflows, and quench the star for- 
mation. This mechanism is also supported by the detection 
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of strong Active Galactic Nuclei (AGN) activity in groups 
and clusters in the form of radio jets (e.g. Boehringer et al. 
1993), or in the distant Universe as seen through quasar 
spectra (e.g. Chartas et al. 2003). 

Semi-analytical models of galaxy formation using cos- 
mological N-body simulations have demonstrated that AGN 
feedback can be responsible for creating a population of red 
and dead massive galaxies (Croton et al. 2006; Bower et al. 
2006; Cattaneo et al. 2006; Somerville et al. 2008). More re- 
alistic cosmological simulations of galaxy formation, includ- 
ing gas dynamics and gas physics, have generally confirmed 
that this mechanism is able to self-regulate the baryon 
content in massive halos to get realistic galaxy masses 
and colours (Sijacki & Springel 2006; Puchwein et al. 
2008; Khalatyan et al. 2008; McCarthy et al. 2010, 2011; 
Dubois et al. 2010, 2011; Teyssier et al. 2011), and that 
it can reproduce the relatively tight Mbh — Mb, Mbh — 
cr b relationships (Sijacki et al. 2007; Di Matteo et al. 2008; 
Booth & Schaye 2009; Dubois et al. 2012). 

Studies of strong gravitational lensing by massive ETGs 
at low redshift (z ~ 0.2) have provided a wealth of in- 
formation about their mass distribution thanks to sizeable 
sample of lenses that have been built up in the last few 
years. The most noticeable result from the SLACS survey 
(e.g. Koopmans et al. 2009; Auger et al. 2010) is that the 
density profile within one half of the projected effective ra- 
dius is found to be close to isothermal (p oc r~ 2 ) with lit- 
tle (~ 10%) scatter around this value (see also Rusin et al. 
2003). It can potentially put some strong constraints on the 
theoretical models as AGN feedback can impact the mass 
distribution of the non-collisional component in the core of 
clusters of galaxies (Dubois et al. 2010, 2011; Teyssier et al. 
2011; Martizzi et al. 2012), even at very early epochs (z ~ 6) 
of the halo formation by blowing away the cold filamen- 
tary gas (Dubois et al. 2013). It is noteworthy to mention 
that these conclusions given by strong lenses apply to low 
(z 0.2) redshift ETGs since most of this work builds 
on the SDSS survey. More recent results have suggested 
that the isothermal profile is also present at higher redshift 
(z ~ 0.5) but natter than at low redshift (Ruff et al. 2011; 
Bolton et al. 2012). 

Deep observations of distant galaxies have re- 
vealed that massive galaxies were more compact in the 
past (Daddi et al. 2005; Trujillo & et al 2006) than they 
are today. These observations challenge the scenario where 
galaxies are built from a monolithic collapse, and favour 
a scenario where mergers drive the formation of mas- 
sive ETGs at z = 0. The size increase of these ob- 
jects due to mergers is expected in semi-analytical the- 
ory and in numerical simulations (Khochfar & Silk 2006b; 
Boylan-Kolchin et al. 2006; Mailer et al. 2006; Naab et al. 
2006, 2007; De Lucia & Blaizot 2007; Bournaud et al. 2007; 
Guo & White 2008; Hopkins et al. 2009; Nipoti et al. 2009; 
Shankar et al. 2013), and is supported by the direct evidence 
of merger remnants of what are the progenitors of today's 
most massive galaxies. Major mergers have been proposed 
as a solution for the evolution of the sizes of galaxies (e.g. 
Khochfar & Silk 2006a), however they are rare events (one 
major merger on average since z — 2) that cannot explain 
the rapid size-mass evolution with redshift (e.g. Bundy et al. 
2009). Dry minor mergers seem to be appropriate candi- 
dates for growing the galaxy sizes as fast as found in the 



observations (Naab et al. 2009; Lackner & Ostriker 2010; 
Oser et al. 2010, 2012; Hilz et al. 2012). 

Hydrodynamical cosmological simulations following the 
formation of elliptical galaxies by Oser et al. (2010), using 
Smoothed Particle Hydrodynamics (SPH), have shown that 
these galaxies can be formed even in the absence of AGN 
feedback, because the fraction of accreted stellar mass dom- 
inates the fraction of in situ-formed stars for low redshift 
massive galaxies. This result has been challenged recently 
by Lackner et al. (2012), using Adaptive Mesh Refinement 
(AMR), where they find that in the absence of AGN feed- 
back, in-situ star formation dominates at all times, even at 
z — 0. The aim of this paper is to probe the role of AGN 
feedback on galaxy morphology, and on the history of mass 
accretion onto galaxies. 

The paper is organised as follow. In section 2, we de- 
scribe our numerical model for galaxy formation and our 
zoom initial conditions. In section 3, we detail the results 
obtained with our simulated galaxies and the particular role 
of AGN feedback on the mass build-up of massive ETGs, 
and we check the consistency of dynamical properties with 
observations. In section 4, we compare our work to previous 
numerical studies solving for gas hydrodynamics. Finally, in 
section 5, we summarise the main results of this work and 
discuss the possible implications. 



2 SIMULATION SET-UP 

2.1 Initial conditions and simulation parameters 

We assume a flat ACDM cosmology with total matter 
(baryons+DM) density Q m = 0.272, baryon density = 
0.045, dark energy density Qa = 0.728, fluctuation am- 
plitude at 8/i _1 .Mpc as = 0.809 and Hubble constant 
H = 70.4km.s _1 .Mpc _1 consistent with WMAP 7-year 
data (Komatsu et al. 2011). 

The simulations are run with the AMR code RAMSES 
(Teyssier 2002). The evolution of the gas is followed using 
a second-order unsplit Godunov scheme for the Euler equa- 
tions. The Riemann solver used to compute the flux at the 
cell interface is the HLLC solver (Toro et al. 1994) and a 
first-order MinMod Total Variation Diminishing scheme is 
applied to reconstruct the interpolated variables from their 
cell-centred values. Collisonless particles (DM, stellar and 
sink particles) are evolved using a particle-mesh solver with 
a Cloud-In-Cell interpolation. 

Two different box sizes Lbox = 100/i _1 .Mpc and 
50/i _1 .Mpc have been used to generate the initial condi- 
tions. Six groups of galaxies in the mass range 4.3 x 10 12 ^ 
Mvir ^ 8.0 x 10 13 M have been selected to be resimu- 
lated at high resolution using the zoom technique, with 
only high resolution DM particles ending up in the virial 
radius of the halos at z = 0. Simulations are allowed to 
refine the initial mesh up to 7 levels of refinement, which 
reach up to a Ax = 1.1 and 0.5 kpc physical length, re- 
spectively, for our low and high resolution initial conditions 
(M D M,re S = 8.2 x 10 7 M and 1.0 x 10 7 M ). A cell is refined 
following a quasi- Lagrangian criterion: if more than 8 dark 
matter particles lie in a cell, or if the baryon mass exceeds 
8 times the initial dark matter mass resolution. 

Six halos are simulated, five out of the 50 /i _1 .Mpc sim- 
ulation box at the highest resolution, and one (the most 
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Table 1. Simulations performed with different sub-grid galactic models and different resolutions, (a) Name of the simulation, (b) Halo 
virial mass at z = 0. (c) Mass resolution of the high resolution DM particles, (d) Stellar mass of the central galaxy at z = 0. (e) Minimum 
cell size, (f) Star formation efficiency, (f) Presence of AGN feedback (Y), or without AGN feedback (N). (h) g— r relative colour magnitude 
of the central ETG accounting for dust extinction, (i) g— r* relative colour magnitude of the central ETG without dust extinction, (j) 
Star Formation Rate (SFR) of the central galaxy at z = 0. 
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massive halo) taken out of the 100 /i _1 .Mpc simulation box 
at lower resolution. These six halos are simulated with the 
same physics (gas cooling, star formation, metal cooling, SN 
feedback, see section 2.2), and with one set without AGN 
feedback (GX simulations), and one set with AGN feedback 
as described in section 2.3 (GXA simulations). Halo masses, 
DM particle mass resolution, and minimum cell sizes are 
summarised in table 1 together with some of the basic prop- 
erties of the central galaxy at z = (mass, colour, star 
formation rate). 

2.2 Physics of galaxy formation 

Gas is allowed to radiate energy away and cool down to 
a minimum temperature of To = 10 4 K due to atomic 
collisions in a gas of primordial composition, i.e. only H 
and He (Sutherland & Dopita 1993). We also account for 
the cooling enhancement due to the presence of metals re- 
leased during Supernovae (SNe) explosions from massive 
stars. Metals are passively advected with the gas, and their 
distribution depends on the history of the gas enrichment 
through SNe explosions. A fixed solar composition of heavy 
elements is assumed for the computation of the gas cooling 
rates. Heating from a UV background is considered follow- 
ing Haardt & Madau (1996) during and after the redshift of 
reionisation z re ion ~ 6. 

Star formation in our model occurs in regions with gas 
number density n > no = 0.1 H.cm -3 using a random Pois- 
son process to spawn star cluster particles, according to a 
Schmidt-Kennicutt law 

where p is the gas mass density, p* the star formation rate 
mass density, e* the star formation efficiency, and tff(p) 
the gas fee- fall time. In these simulations, we set the effi- 
ciency of star formation to a low value e* = 0.02 that is 
in good agreement with observational surface density re- 
lationships of galaxies (Kennicutt 1998), and local giant 
molecular clouds (Krumholz & Tan 2007). Each star clus- 
ter particle has a mass resolution of ra* = poAx 3 (po is 



the gas mass density threshold of star formation) where a 
random Poisson process is used to create a new star parti- 
cle (Rasera & Teyssier 2006; Dubois & Teyssier 2008). The 
stellar mass resolution reaches 5.1 x 10 5 Mq for our most 
resolved simulation with Ax = 0.5 kpc. We ensure that no 
more than 90% of the gas in a cell is depleted during the 
star formation process for numerical stability. 

We account for the mass and energy release from type II 
SNe assuming a Salpeter Initial Mass Function (IMF). Using 
this IMF, 10% of the massive stars above 8Mq end their 
life into a type II SN releasing 10 51 erg per 10 M©. For each 
individual SN explosion, after 10 Myr, energy is released 
spherically in a kinetic mode transferring mass, momentum 
and energy with a mass loading parameter f w = 1 that 
mimic a Sedov blast wave explosion (see Dubois & Teyssier 
2008, for further details). We assume that type II SNe release 
all their mass into the gas (no stellar remnant forms) with a 
y = 0.1 stellar yield, which is the fraction of primordial gas 
transformed into heavy elements and released back into the 
gas. With this prescription, we do not take into account the 
energy and mass release from stellar winds ( AGB stars) , nor 
from long-lived type la SNe. 

The gas follows an adiabatic equation of state with an 
adiabatic index 7 a diab = 5/3. In order to take into account 
the thermal impact of the heating of the gas by SNe, we 
modify the temperature at high gas density p > po with a 
polytropic equation of state 

T = To(£f~\ (2) 

where p is the polytropic index of the gas. The adopted 
value of p — 4/3 is comparable to the value obtained 
in Springel & Hernquist (2003) considering the multiphase 
structure of the ISM with stellar heating. We point out 
that this value of p = 4/3 does not rigorously ensure that 
gas fragmentation is avoided due to numerical instabilities 
(Truelove et al. 1997), as the Jeans length is proportional to 
the gas density 

Aj = 10 -KoTH^)" 1/3kpC - (3) 
This last formula shows that at very high gas densities, the 
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Jeans length can be s" mailer than our minimum resolution 
and would provide spurious fragmentation of the gas. How- 
ever, the gas cannot indefinitely condense because of the 
finite force resolution sampling and the star formation pro- 
cess that removes gas in high density regions. 



2.3 AGN feedback 

Our modelling of AGN feedback closely follows the im- 
plementation detailed in Dubois et al. (2012). We re- 
call here the main aspects of this modelling, and refer 
to Dubois et al. (2010) and Dubois et al. (2012) for the tech- 
nical issues. Black hole particles are modelled with the sink 
particles technique from Bate et al. (1995), and particu- 
larly Krumholz et al. (2004) for mesh codes. The BHs are 
created in the centres of galaxies, with a sufficiently large 
exclusion radius (100 kpc) to prevent formation of multiple 
massive BHs per galaxy. The initial BH seed mass is cho- 
sen to be 10 5 M . BHs can merge when they are sufficiently 
close to one another (< 4 Ax), determined using a Friend- 
Of-Friends algorithm on sink particles performed on-the-fly. 

BHs can grow by accretion at a Bondi-Hoyle-Lyttleton 
rate (Bondi 1952) 



M BH 



47raG 2 M^ u p 
(cf+n 2 ) 3 / 2 



(4) 



where G is the gravitational constant, p the average gas 
mass density, c s the average sound speed, u the average gas 
velocity relative to the BH, and a a dimensionless boost 
factor with a = (p/po) 2 when p > po and a = 1 other- 
wise (Booth & Schaye 2009) in order to account for our in- 
ability to capture the cold and high density regions of the 
gas at these kpc resolutions. 

The accretion rate onto BHs is limited by its Eddington 

rate 



M Ed d 
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€ r (JTC 



(5) 



where <7t is the Thompson cross-section, c is the speed of 
light, m p is the proton mass, and e r is the radiative efficiency, 
assumed to be equal to 0.1 for the Shakura & Sunyaev 
(1973) accretion onto a Schwarzschild BH. 

A drag force is applied to the BH particles in order 
to replicate the action of dynamical friction of the gas on 
massive BHs that happens on unresolved parsec scales. The 
dynamical friction force is defined as 



p f a -f GM BH \ 2 
Fof = / gas 47rap I — I . 



(6) 



where / gas is a fudge factor whose value is between and 
2 and is a function of the Mach number Ai = u/c s < 
1 (Ostriker 1999; Chapon et al. 2011), and where a is the 
boost factor already introduced in equation (4). This ex- 
tra force term for BH dynamics has been successfully intro- 
duced in high redshift simulations with ~ 10 pc resolution to 
maintain BHs close to the centre of the galaxy (Dubois et al. 
2013). 

In this paper, we assume that AGN feedback proceeds 
uniquely in the jet mode powering large bipolar outflows. 
The total energy released by the AGN is 
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Figure 2. Stellar mass of central massive galaxies as a function 
of their halo mass at z = for the simulations without AGN 
(black), and for the simulations with AGN feedback (red). The 
solid curves are the observational fit from Moster et al. (2010) ei- 
ther assuming a Kroupa (black) or a Salpeter (blue) IMF, with 
the 1, 2 and 3 a standard deviation (gray shaded areas). Cyan 
solid lines indicate the amount of stars formed at a constant effi- 
ciency / C onv = M s /f] D M h . The presence of AGN feedback reduces 
the total amount of stars in the central massive galaxy by a factor 
7. 



where ef = 1 is a free parameter chosen to reproduce the 
MeH-Mb, MBH-o"b, and BH mass density in our local Uni- 
verse (see Dubois et al. 2012). For this jet mode of AGN 
feedback, mass, momentum and kinetic energy are continu- 
ously released in the surrounding gas along a jet axis defined 
by the local angular momentum of the gas, and such that the 
jet velocity is 10 4 km.s _1 , or equivalently the mass loading 
factor of the jet is 77 = 100. 



3 RESULTS 

3.1 Changing the baryon content and the 
morphologies of massive galaxies 

The basic picture of the formation of massive halos is that 
they form at high redshift by cold filamentary infall driv- 
ing high levels of star formation, and at low redshift by a 
fainter diffuse (and hot) accretion, as illustrated for the halo 
G5A in Fig. 1 (Birnboim & Dekel 2003; Keres et al. 2005; 
Ocvirk et al. 2008; Dekel et al. 2009). For groups of galax- 
ies, satellites represent a non-negligible fraction of the total 
stellar mass in the halo with the large number of satellites 
orbiting inside the virial radius of the halo. Note that some 
of these satellites have no dense star-forming gas compo- 
nent. This gas has been stripped away by the ram pressure 
of the hot diffuse gas as can be seen for some satellites with 
tails of gas. 

Massive halos usually suffer from over-cooling of their 
gas content in absence of any strong feedback prescription. 
Vanilla hydrodynamical simulations (i.e. without any strong 
feedback) produce large amount of stars in massive central 
galaxies as shown in Fig. 2. The measured stellar mass M s 
is given by the total stellar mass of the central galaxy as 
detected by the AdaptaHOP algorithm (Aubert et al. 2004; 
Tweed et al. 2009) performed on stars, and, thus, captures 
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Figure 1. Top and middle rows: Images of the projected gas number density at z = 5 (top left), z = 3 (top centre), z = 1.5 (top right), 
z = 1 (middle left), z = 0.5 (middle centre), and z = (middle right) around the most massive galaxy of the G5A simulation; and 
bottom: Images of the projected DM surface density (left), and stellar surface density (right) at z = 0. The large circles correspond to 
the virial radius of the host halo and the small circles to 10 per cent of the virial radius. The image sizes are in physical length. The gas 
distribution around the central galaxy shows a complex filamentary network at high redshift, and a smooth spherical distribution of gas 
with a large number of satellites at intermediate and low redshift. 
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Figure 3. Average conversion efficiency of the central galaxies 
/conv = M s //bMh as a function of redshift for the simulations 
including AGN feedback (red) or without AGN feedback (black) 
with the standard deviation (error bars). The values are averaged 
over the main progenitor of the central galaxy for the six halos. 
The conversion efficiency of stars in the central galaxy is strongly 
reduced by the presence of AGN feedback and the reduction takes 
place already at high redshift (a factor 2 reduction at z = 4). 



the stellar content in the core of the halo together with its 
related intra-cluster light, and removes all satellites (consid- 
ered as "sub-galaxies" of the central most massive galaxy). 
The halo mass Mh is the total mass of gas, stars, DM (and 
BHs if any) enclosed within the virial radius provided by 
the halo finder algorithm (using AdaptaHOP on DM parti- 
cles). In the absence of AGN feedback, stars are converted 
in the central galaxy with a / con v = 0.7 global efficiency, 
where / con v = M s //bMh and fh = Qb/Qm = 0.165, while 
observations (e.g. from Moster et al. 2010) suggest that the 
value for groups of galaxies is of 0.1 (for a halo mass of 
Mh ~ 10 13 M©) and decreases with increasing halo mass. 
AGN feedback manages to reduce the stellar mass of the 
central massive galaxies at z — by a factor 7 compared to 
the no AGN feedback case, but these galaxies still remain 
too massive, in particular even for those belonging to the 
most massive halos. 

Note that if, as suggested by recent observa- 
tions using lensing (or dynamics) and photome- 
try (Grillo et al. 2008; Grillo & Gobat 2010; Treu et al. 
2010; Auger et al. 2010; Spiniello et al. 2011; Thomas et al. 
2011; Sonnenfeld et al. 2012; Cappellari et al. 2012) or 
spectroscopy (Spiniello et al. 2012; Conroy & van Dokkum 
2012), the IMF is considered to be more bottom- heavy 
for the most massive galaxies, then the observational 
relation of M s versus Mh should be pushed towards higher 
values of M s . Assuming a Salpeter (1955) IMF instead of 
a Kroupa (2001) (or Chabrier 2003) IMF, the ratio M s /M h 
increases by +0.3 dex, which puts our values of the central 
stellar mass for the four most massive galaxies (with AGN 
feedback) in better agreement with observations, but in 
less better agreement for the two least massive galaxies 
(G1A and G2A). It exists a possible smooth transition 
from a Kroupa to a Salpeter IMF for galaxies with velocity 
dispersion around 100 — 200km.s _1 (Cappellari et al. 
2012), values which correspond to the typical ones of the 



galaxies at z = when simulated with AGN feedback (see 
section 3.3). 

Fig. 3 shows the evolution of the stellar conversion effi- 
ciency /conv averaged over the six central galaxies as a func- 
tion of redshift. In the absence of AGN feedback, / con v is con- 
tinuously increasing with decreasing redshift. Large amounts 
of gas are piling up in the central galaxy due to over-cooling, 
but, as it takes several free-fall times (tff/e* = 50£ff ~ 1 Gyr 
for a typical gas number density of n — 10 cm -3 in star 
forming regions, see Fig. 1) to convert the cold gas content 
into stars, the value of / con v is lower at high redshift. The 
same trend for the evolution of /conv is observed at high red- 
shift above z ~ 2 for the simulations with AGN feedback. 
However, during the early phase of galaxy formation, the 
presence of AGN feedback allows for a rapid quenching of 
the star formation with a reduced value of / con v by a factor 
2 at z — 3 — 4 compared to galaxies simulated without AGN. 
From z ~ 2 to z = 0, the value of / con v slowly decreases by 
a factor 2. While halos continue to grow by diffuse accretion 
and mergers, the gas is efficiently prevented to cool further 
down by AGN feedback, and the central galaxy stops ac- 
tively forming stars (see the values of the SFR obtained at 
z = in table 1). 

By simulating the stellar light emission (using the SDSS 
u, g and r filter bands) of the central galaxies, we observe 
that, because of AGN feedback, they become less luminous, 
more red and with rounder shapes (see Fig. 4). In the ab- 
sence of AGN feedback, the central galaxies exhibit a mas- 
sive blue disc component in their centres, that is confirmed 
by the measured colours, around g — r — 0.45 without 
dust extinction. For these over-massive objects, there is a 
strong dust lane component (not represented in the images 
of Fig. 4), which reddens them by A(g — r)d us t = 0.2 (see 
table 1). In contrast, simulated galaxies with AGN feedback 
have colours between 0.60 < g — r < 0.73 (no dust extinc- 
tion), and for some of them a contribution is adopted for 
dust extinction of A(g — r)dust = 0.1 (closer to zero for some 
others) . We get better consistency of galaxy colours with ob- 
servations when AGN feedback is considered (Bernardi et al. 
2005). 

We measure the average mass fractions of DM, stars and 
gas within the 0.1r v ir of the halo as a function of redshift 
for the simulations with or without AGN feedback (Fig. 5). 
We find that the fraction of DM mass is more important in 
the core of halos when AGN feedback is present, going from 
(/dm) = 0.53 without AGN up to (/dm) = 0.83 with AGN 
at z = 0. The DM fraction is constant over time in the ab- 
sence of AGN feedback, but increases with time for the AGN 
case. Note that the DM fraction is already higher in the AGN 
case at high redshift compared to the no AGN case, but this 
difference is enhanced at low redshift. Without AGN feed- 
back the fraction of stars is increasing with time and the gas 
fraction is decreasing, while the baryonic fraction is constant 
(DM fraction is constant). Thus, galaxies are gas-rich (wet) 
at high redshift, and gas-poor at low redshift because they 
consume the gas efficiently in the star formation process. 
Both the stellar and the gas mass fractions are reduced due 
to the presence of AGN feedback, from (/ s ) = 0.42 down 
to (f B ) = 0.15, and from (/ g ) = 0.04 down to (/ g ) = 0.016 
at z = 0. The difference in stellar and gas mass fractions is 
also seen at higher redshift, but in smaller proportions than 
at low redshift. Finally, AGN feedback already impacts the 
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Figure 4. Stellar emission of the central galaxies from Gl to G6 (left to right) without (top) or with AGN feedback (bottom) at z = 0, 
as they would be observed in u-g-r filter bands. Extinction by dust is not taken into account in these images. Arbitrary units are used 
for stellar emission but with similar minima and maxima for all the first fifth columns, and with decreased intensities by 0.5 dex for 
the last column to avoid image saturation. The size of each panel is 140 kpc. Galaxies without AGN feedback exhibit a massive blue 
disc component with a halo of stars extending to large distance. AGN feedback reduces the total amount of emitted light, and galaxies 
appear more early type with a weak disc component for some of them. 



mass content in galaxies early on (at z ~ 4), but its effect is 
largest at z = 0. 

Thanks to an exquisite resolution, we are able to go be- 
yond the global properties listed above and investigate the 
internal dynamical structure of galaxies. The morphology of 
galaxies is also transformed by the effect of AGN feedback 
as shown in Fig. 6. Circular velocities v c — v / CM to t(< r)/r 
exhibit a peak close to the centre of the galaxy, below 5 kpc, 
when AGN feedback is not included. This peak is produced 
by the large concentration of material in the central region of 
galaxies, where the total matter distribution is largely dom- 
inated by stars (see Fig. 7), and is a characteristic feature of 
the overcooling problem of baryons (e.g. Scannapieco et al. 
2012; Few et al. 2012). In the opposite case, for the sim- 
ulations including AGN feedback, the peak in the circular 
velocity profiles disappears and profiles become flat (isother- 
mal) up to very large distances away from the centre, with 
values close to the circular velocity measured far away from 
the centre (at rsoo)- With the presence of AGN feedback, 
the maximum of the circular velocity curves decreases, the 
contribution from stars to circular velocity is strongly at- 
tenuated, and replaced by its DM component (Fig. 7). Note 
that the contribution from gas to the circular velocity curves 
at z = is negligible with or without AGN feedback, which 
suggests that these massive galaxies are relatively gas-poor 
independently of the presence of AGN. We notice that the 
ratio of radial velocity dispersion over the circular velocity 
is roughly constant and close to a factor y/2, therefore, it 
gives further evidence for an isothermal profile. It assumes 
that orbits are isotropic, i.e. the velocity tensor is close to 
isotropic. Further investigations along this line are left for 
future work. 

AGN have also important consequences for the total 
amount of rotation and velocity dispersion in these massive 
galaxies. Without AGN feedback, the stellar component of 
the massive galaxies at z = is supported by the rotation of 
stars, with a comparable but lower contribution from disper- 
sion, particularly within the half mass radius of the galaxy. 
Such galaxies would be classified as fast rotators as the v/a 
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Figure 5. Average fraction of DM (solid), stellar (dotted), and 
gas (dashed) mass within 0.1r v j r as a function of redshift for the 
six halos with (red) or without (black) AGN feedback. The error 
bars are the standard deviation. The fraction of DM mass is in- 
creased, and the fractions of gas and stellar mass are decreased 
by the effect of AGN feedback. 

ratio for stars approaches one (even though their masses and 
colours are unrealistic). At a close distance from the centre, 
the velocity support is dominated by the dispersion (com- 
pact bulge of stars), and also far away, in the intra cluster 
light. The picture is changed when AGN feedback is present: 
velocity dispersion dominates at any distance from the cen- 
tre at z = 0. The disc component has considerably shrank 
and, now, is barely detectable below 5 kph at z = (Fig. 7). 

Fig. 8 shows the average v/a ratio, where v is the mass- 
weighted rotational velocity of stars and a is the mass- 
weighted velocity dispersion of stars measured within r e ff, 
where r e ff is the effective radius at which half of the pro- 
jected stellar mass is enclosed. Both quantities are integrated 
quantities of the projected distribution of stars, and the v/a 
ratios showed in Fig. 8 are averaged over the six galaxies 
(with or without AGN feedback). Without AGN feedback, 
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Figure 6. Circular velocities averaged over the six halos with 
(red) or without (black) AGN feedback as a function of radius 
at z = 0. Velocities are renormalised by the circular velocities at 
rsoo- Shaded areas correspond to the standard deviation. Simula- 
tions without AGN feedback exhibit a central peak in the circular 
velocity that is characteristic of too large mass concentration in 
the centre of the galaxy, while simulations with AGN feedback 
show flat velocity curves with a value close to the circular veloc- 
ity at rsoo- 

the galaxies have large values of v/a above 1, i.e. are dom- 
inated by the rotation of stars at all redshifts. Due to the 
increased proportion of accreted material (see next section) , 
the v/a ratio decreases with time and tends towards 1 at 
z = 0. Thus the massive galaxies at z = simulated with- 
out AGN feedback are disc-like galaxies. The transformation 
of rotationally-supported discs into dispersion-dominated el- 
lipsoids by the presence of AGN feedback is observed for all 
our simulated central galaxies at low redshift (z < 1). The 
same evolution is found for galaxies simulated with AGN 
feedback as without AGN feedback: the ratio of v/a also 
decreases with time, but AGN transforms the galaxies into 
slow rotators (v/a < 1) where galaxies become essentially 
supported by the velocity dispersion of the stars, with a 
value slightly above 1 at z = 4 down to a value of 0.5 at 
z = 0. The importance of rotation over dispersion is already 
reduced by a factor 2 at high redshift (z = 3 — 4). 

We interpret this morphological change as the signature 
that massive galaxies are transformed from systems domi- 
nated by in situ star formation (without AGN feedback) 
into systems dominated by accretion of satellites (with the 
presence of AGN). We will now clearly demonstrate this 
mechanism in the following section. 

3.2 From in situ to accretion-dominated systems 

The fraction of stars locked into satellites is modified by the 
presence of AGN as seen in Fig. 9. Without AGN feedback 
the fraction of stars in satellites M s , S at/M s ,tot is increasing 
with time (M SjSa t is the mass of stars in satellites within 
r v ir, and M s ,tot is the total mass of stars within r v ir), and 
reaches values of 20-30 per cent at z — 0. It shows that satel- 
lite galaxies are more present in the local Universe than in 
the early Universe, but, in the absence of AGN feedback, the 
total stellar mass in the halo is still dominated by the central 
galaxy. When AGN feedback is active, it efficiently reduces 



1000 




(kpc) 



Figure 7. Illustration of the morphological transformation of the 
central massive galaxy of the halo of G3 at z = due to AGN 
feedback. Both panels show the circular velocities as a function 
of radius (black) of all components (solid), stars (dashed), DM 
(dotted), gas (dot-dashed), and BHs (triple dot-dashed) without 
AGN feedback (top panel), or including AGN feedback (bottom 
panel). The measured rotational velocities for stars ^ s ,rot (blue 
plus) and radial velocity dispersions c srad (orange triangle) are 
also indicated. The 3D stellar half-mass radius 7*h a if an d 10 per 
cent of the virial radius 0.1r v j r are indicated as vertical bars. The 
circular velocity is decreased by a factor 2 due to AGN feedback. 
The circular velocity is dominated by the stellar component up to 
large distance (30 kpc) in the AGN case, while it dominates only 
in the very central region (7 kpc) in the AGN case, and the DM is 
the principal contributor outside these characteristic radius. The 
galaxy is transformed from a rotating disc (no AGN case) into a 
velocity dispersion- dominated ellipsoid (AGN case). 



the total stellar mass in the central galaxy (Fig. 3) by blow- 
ing gas away and preventing accretion. In consequence, the 
contribution of satellites to the total stellar mass within the 
virial radius of the halos increases such that the satellites 
represent half of the stellar mass at z = 0. The increased 
fraction of mass locked into galaxy satellites due to AGN is 
a clear indicator that the mass build-up of central galaxies 
should have an increased proportion of stellar mass gained 
through accretion. 

To demonstrate this effect, we define the in situ formed 
stars as the stars formed within the main progenitor of the 
galaxy (within 10 per cent of the virial radius of the host 
halo), and the accreted stars are defined to be the stars 
formed outside of the main progenitor of the galaxy and 
acquired through capture of satellites or by diffuse accre- 
tion. Fig. 10 shows the amount of stellar mass in the cen- 
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Figure 8. Average v/cr ratio as a function of redshift for the 
central galaxy of the six simulated halos with (red) or without 
(black) AGN feedback, with the standard deviation (error bars). 
The galaxies simulated without AGN feedback are fast rotators 
at any redshift, while galaxies simulated with AGN feedback are 
slow rotators. 
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Figure 9. Average fraction of the total stellar mass in the halo 
locked into satellite galaxies without AGN feedback (black) and 
with AGN feedback (red). The values are averaged over the six 
simulated halos, and the error bars correspond to the standard 
dispersion. With the presence of AGN feedback, the fraction of 
stars locked into the galaxy satellites is increased with almost half 
of it in the satellites at z = 0. 



tral massive galaxy formed in situ and the amount of stars 
that are accreted, for the different simulations, depending 
on the presence of AGN feedback. Galaxies are dominated 
by the in situ growth of the stellar mass at high redshift 
due to the large fraction of baryons available in the form 
of cold star-forming gas (Tacconi et al. 2010; Daddi et al. 
2010; Dubois et al. 2012). As a consequence, the specific star 
formation rates (star formation rate per galaxy stellar mass) 
are maintained at high values, even when AGN feedback 
is included (Kimm et al. 2012). Conselice et al. (2012) have 
shown that observed massive galaxies M s > 10 11 M© acquire 
the majority of their stars 61 ± 21 % through the in situ star 
formation between redshift 1.5 < z < 3 which is consistent 
with our findings. As time goes on, the gas reservoir is con- 



sumed in the star formation process and ejected through the 
galactic winds produced by AGN feedback, and the specific 
star formation rates decline. The result is that mergers have 
an increasing contribution to the mass build-up of these mas- 
sive galaxies (De Lucia & Blaizot 2007), and the fraction of 
in situ-formed stars declines with time (decreasing redshift) 
as seen for all the six simulated halos independently of the 
presence of feedback. 

The difference between the two simulation types (with 
or without AGN) is that, for the "vanilla" model (no AGN), 
the central galaxy mass is dominated by the in situ star for- 
mation (stars are formed locally) at z — with /msitu — 55 
per cent of the total stellar mass on average, while, for the 
AGN feedback runs, the stellar mass is dominated by the 
accretion of satellites at z — with /msitu = 40 per cent on 
average, as seen in Fig. 11. It is important to stress that both 
the mass of in situ formed stars and the mass of accreted 
stars at z = are decreased because of the presence of AGN 
feedback (Fig. 10), but AGN feedback quenches the star for- 
mation more efficiently in the central galaxy than it does in 
the galaxy satellites. Note that, as suggested by Fig. 3, the 
total mass of stars in the central galaxy is already affected 
by the presence of AGN feedback early on {z = 3 — 4), and 
that the fraction of in situ-formed stars is decreased early 
on when the activity of AGN is the strongest (Dubois et al. 
2012). 

As seen in the previous section (see Fig. 7 and Fig. 8), 
the galaxies are transformed from rotationally supported 
(discs) into dispersion dominated (ellipticals) systems. We 
will show that this is a direct consequence of the trans- 
formation of galaxies being build-up by in situ star forma- 
tion (without AGN) into accretion dominated systems (with 
AGN), and that the transformation has important implica- 
tions for the dynamical properties of these objects. As an 
illustration, Fig. 12 shows the stellar density profiles of the 
central galaxy G4 and the central galaxy G4A at z = 0, 
decomposed into in situ and accretion components. The in 
situ component dominates the total stellar density profile in 
the core of the galaxy, whereas the accreted material dom- 
inates in the outskirts. This is a natural outcome of the in 
situ stars forming in the depth of gravitational potentials, 
and the satellites being tidally stripped while falling into 
halos on elliptical orbits, thus depositing accreted stars far 
away from the centre (Naab et al. 2009; Peirani et al. 2010; 
Hilz et al. 2012). 

AGN feedback changes the radius at which the transi- 
tion between an in situ-dominated density towards an accre- 
tion dominated density occurs. As the proportion of accreted 
material is increased due to the efficient quenching of the in 
situ star formation by AGN feedback, the radius at which 
accreted material dominates is pushed towards the interior 
part of the halo (e.g. from 24 kpc to 14 kpc for the G4X 
halo at z = 0, see Fig. 12). This is an important concept: as 
the effective radius (or half-mass radius for 3- dimensional 
profiles) associated with the accreted material is larger than 
those of the in situ component, the result is that galaxies get 
puffed up due to the presence of AGN feedback. 

Fig. 13 shows the evolution in redshift of the effective 
radius, and the in situ and accreted effective radii for all 
halos with or without AGN feedback. All simulations show 
a larger value of their effective radii in the cases with AGN 
feedback, which is a direct consequence of the mechanism 




Figure 10. Top: Fraction of in situ formed stars and, bottom: total stellar mass (solid lines), mass of in situ formed stars (dotted lines) 
and of accreted stars (dashed lines) as a function of redshift for Gl (top left), G2 (top middle), G3 (top right), G4 (bottom left), G5 
(bottom middle), and G6 (bottom right). The results for simulations without AGN feedback are in black (GX) and with AGN feedback 
in red (GXA). The grey horizontal lines indicate the 50 per cent level of in situ versus accreted stellar mass. 



described above: as accretion has a more important role, 
galaxy sizes are increased. This mechanism is also confirmed 
by the fact that high redshift effective radii are closer to 
their in situ effective radii, while low redshift effective radii 
are closer to the accreted effective radii (respectively when 
in situ, or accreted mass dominates, see Fig. 10). The aver- 
age logarithmic slope of the evolution of the effective radius 
with redshift d log r e s / dz measured between redshift and 
2 is equal to —0.23 in the case without AGN feedback, and 
to —0.39 in the case with AGN feedback. Thus, the effec- 
tive radius of galaxies simulated with AGN feedback evolves 
faster with time than galaxies without AGN. This is a direct 
consequence of the galaxies having an increased fraction of 
accreted material in their mass assembly in the former case. 

The effective radius of the accreted material is increased 
for all galaxies with the presence of AGN feedback between 
^ z < 1 (dashed lines in Fig. 13), where the sharp tran- 
sition from in situ to accreted dominated galaxies occurs 
(Fig. 11). The explanation is that tidal stripping becomes 
more efficient at removing stars from satellites. The charac- 
teristic radius above which tidal stripping can remove mate- 



rial from a satellite is (King 1962) r tB = £>(0.5M sa t /M s ) - 5 
where D is the distance between the central galaxy and the 
satellite. The distance at which satellites start losing stars is 
when rts < Teff, or D < r e ff (2M s /M S at) " 5 . Thus, if galaxies 
are less compact (i.e. they have larger effective radii), satel- 
lites will be more subject to tidal stripping early on, and will 
lose their stars at larger distances. This is compatible with 
the fact that galaxies are less compact with AGN feedback 
as observed from our simulations. 

The behaviour of the effective radius of the in situ com- 
ponent together with the effect of AGN feedback shows more 
variations (dotted lines in Fig. 13). Some simulated galax- 
ies show an increase of their in situ effective radius, while 
some other have a more compact in situ stellar content. 
We identify two effects that can act in opposite directions. 
As galaxies build from inside out in the classical picture of 
gas accretion, galactic discs of gas tend to be smaller due 
to strong AGN feedback as the accretion is quenched onto 
galaxies, an effect which can potentially decrease the size 
of the in situ component. On the other hand, strong bursts 
of feedback have been identified as creating strong varia- 
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Figure 11. Average fraction of in situ formed stars as a function 
of redshift without AGN feedback (black) and with AGN feedback 
(red). The grey horizontal lines indicate the 50 per cent level 
of in situ versus accreted stellar mass. The error bars indicate 
the standard deviation to the average value. The in situ fraction 
of stars decline with time, and the presence of AGN feedback 
reduces this fraction further down to turn galaxies into accretion 
dominated systems at z = 0. 
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Figure 12. Total stellar density profiles (solid lines) as a function 
of radius at z = for the central galaxy of G4 (black) and G4A 
(red), with the in situ selected stars (dotted lines) and accreted 
stars (dashed lines). The vertical bars indicate the position of 
the radius of half mass for the different components. The half- 
mass radius of the accreted component is larger than the half- 
mass radius of the in situ component. As the accreted material 
dominates the stellar density profiles at smaller radii when AGN 
feedback is turned on, the half-mass radius of the total stellar 
mass (in situ + accreted) increases. 

tions in the gravitational potential well, and can smooth 
the density profiles of both stars and DM in the most central 
parts (Peirani et al. 2008). This effect can increase the size 
of galaxies and DM halo cores, but depends on the details of 
the accretion history onto the central galaxy (Martizzi et al. 
2012). Some of the galaxies show a larger in situ effective ra- 
dius in the AGN case compared to the no AGN case at high 
redshift where the wet and chaotic accretion takes place. 
This can trigger strong bursts of AGN feedback that may 
be effective at flattening the density profiles by star forma- 



tion quenching. Also, we observe such flattening in the total 
mass density profiles (and in the DM component) and we 
will discuss this further in the next section by comparing 
our results from simulations to observations. 

3.3 Implications for galaxy scaling relations and 
mass distribution 

The scaling relations between mass and kinematics can put 
strong constraints on the modelling of the population of mas- 
sive galaxies. In Fig. 14, we show the relations between the 
effective radius r e ff and the stellar mass of the central galaxy 
at different redshifts. Independently of the adopted IMF, 
all galaxies simulated without AGN feedback lie below the 
r e ff-M s observational relation at z = 0.2 from the SLACS 
data (Auger et al. 2010) that consists of ETGs within the 
mass range 10 11 < M s < 6 x 10 11 M . These galaxies, suf- 
fering from over-cooling and being too massive, are also 
too compact. This increased compactness is expected if the 
amount of in situ star formation is too large, as the in situ 
formed stars have small effective radius compared to the 
effective radius of the accreted stellar distribution (Fig. 13). 

The presence of AGN feedback increases the effective 
radius, while at the same time it reduces stellar masses 
and greatly improves the agreement between the simula- 
tions and observations for the four most massive galaxies 
(G3A, G4A, G5A, G6A). The simulation points with AGN 
feedback stand close to the observational result for the four 
most massive galaxies at z — 0.2, within the 2a error bars 
(for a Chabrier IMF). If a Salpeter IMF is assumed, the 
AGN size-mass points of the four most massive galaxies at 
z = 0.2 stand within the 5cr relation. The two least massive 
galaxies are largely above the observational relation, which 
makes these galaxies too puffy. Such galaxies exist but are 
rare objects in the real Universe (Williams et al. 2010). It is 
difficult to explain this particular behaviour, but one should 
note as they are the lightest galaxies, they are also the galax- 
ies that are the most affected by the finite resolution. At 
redshift z = 2, the effective radius of these galaxies is very 
close to the resolution limit of the grid (r e s — 2 kpc against 
a Ax — 0.5 kpc resolution), and the galaxy gets extra sup- 
port against gravity from the finite cell size. At even higher 
redshifts, this effect is true for all simulated galaxies as they 
are smaller, and could have some non-negligible effects on 
the final size-mass relation at z = 0. 

Galaxies are evolving in size with time while they accu- 
mulate stars (Fig. 13), and there is a small variation of the 
size- mass relation with redshift for the AGN feedback sim- 
ulations, with galaxies being more compact at high redshift 
(i.e. the distribution of r e ff versus M s at high redshift is be- 
low this relation at low redshift, Fig. 14). Observations sug- 
gest such an evolution with time in the compactness of galax- 
ies (Trujillo et al. 2007; Cimatti et al. 2008; Saracco et al. 
2009; van Dokkum et al. 2010; Huert as- Company et al. 
2013), and the effective radius is of the order of ~ 1.5 — 2 
smaller at redshift 1 than at redshift at constant mass 
for ETGs (Huertas- Company et al. 2013). On the opposite, 
simulations points without AGN feedback show the oppo- 
site trend: galaxies tend to be less compact at high redshift. 
Note that there is likely an effect of limited resolution on 
the size evolution, that more importantly impacts the low- 
mass galaxies in their early stages (see the flattening of the 




Figure 13. Effective radius of the central galaxy (solid), of the in situ component only (dotted) and of the accreted component only 
(dashed) as a function of redshift for the simulations GX without AGN feedback (black) or GXA with AGN feedback (red). All galaxies 
included AGN feedback show an increased value of the effective radius relative to the no AGN case. 




Figure 14. Effective radius r e ff versus stellar mass M s for our 
six halos Gl (crosses), G2 (stars), G3 (pluses), G4 (squares), G5 
(diamonds), and G6 (triangles) with (red) or without (blue) AGN 
feedback computed at different redshift s: z = 0.2 are large sym- 
bols and z = 1 are small symbols. The observed relations from the 
SLACS (at z ~ 0.2) are overplotted as solid lines with their 2a 
standard deviation (dot) assuming either a Chabrier (black) , or a 
Salpeter IMF (green). Galaxies simulated with AGN feedback are 
less compact than those simulated without AGN feedback, and 
are more compact at high redshift than they are at low redshift. 



curves at z > 1 in Fig. 14). Because of additional numerical 
support we can expect their size to be artificially increased 
when they are close to the resolution limit. 

We also measured the projected velocity dispersion of 
stars a s within r e ff/2 relative to their galaxy mass M s 
(Fig. 15). The galaxies without AGN feedback show large 
values of a s . This is the direct consequence of the size-mass 
relation: more compact galaxies need more dynamical sup- 



Figure 15. Same as Fig. 14 for the velocity dispersion of stars 
a s within r e ^/2 versus stellar mass M s . Galaxies simulated with 
AGN feedback have lower velocity dispersion values than for 
galaxies simulated without AGN feedback. 

port (velocity dispersion). At z = 0.2 the four most massive 
galaxies without AGN lie above 2a of the observational re- 
lation, while the two least massive lie within the relation 
(for a Chabrier IMF). As can be predicted from the size- 
mass relation, the presence of AGN feedback produces more 
extended galaxies and, thus, galaxies with lower values of 
stellar velocity dispersion. The four most massive galaxies 
with AGN have values of a s -M s in better agreement with 
observations than in the absence of AGN. As for the size- 
mass relation, the a s -M s points for the two least massive 
galaxies are at more than 2a away from the observational 
relation at z = 0.2 if AGN feedback is incorporated. These 
low mass ETGs are above the size-mass observational fit and 
since they need less kinematic support because they are too 
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Figure 16. Stellar mass M s -plane at z = 0.2 for galaxies with 
AGN feedback (red) and galaxies without AGN feedback (black). 
The observed relations from the SLACS (at z ~ 0.2) are overplot- 
ted as solid lines with their 2a standard deviation (dot) assum- 
ing either a Chabrier (black), or a Salpeter IMF (blue). Galaxies 
simulated without AGN feedback lie above the M s -plane, while 
galaxies with AGN feedback are closer to M s -plane (favourising 
a Chabrier IMF). 
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Figure 17. Average slope of the total mass density profile as a 
function of redshift for the six halos with (red) or without (black) 
AGN feedback. The slope is measured between r e ^/2 and r e ^ 
and on the average mass density within radius r. Error bars cor- 
respond to the standard deviation. Observations from the SLACS 
sample (Auger et al. 2010) at z = 0.2 are plotted in cyan. The 
observational trend of the slope with redshift found by Ruff et al. 
(2011) and Bolton et al. (2012) are plotted in violet and orange 
respectively with their intrinsic scatter. AGN feedback flattens 
the total density profiles. AGN feedback produces shallower total 
density profiles within galaxies. 

extended, they lie below the velocity dispersion-mass obser- 
vational law. 

An even more stringent observational constraint to 
reproduce is the fundamental plane of ETGs. Here, 
we measure the mass equivalent, namely the Ms- 
plane (Bolton et al. 2007, 2008; Auger et al. 2010), that 
links the galaxy size to the galaxy mass and velocity dis- 
persion (Fig. 16). All galaxies simulated without AGN feed- 
back stand above the M s -plane. For the AGN feedback sim- 



Figure 18. Same as Fig. 14 for the DM fraction /dm within r e ff/2 
versus the velocity dispersion of stars a s . Galaxies simulated with 
AGN feedback have larger DM fraction within r e ^/2 than galaxies 
simulated without AGN feedback. 

ulations, we find that the four most massive galaxies stand 
closer to the M s -plane than in the case without AGN feed- 
back and they favour a Chabrier IMF. For the two least 
massive galaxies including AGN, they stand below the M s - 
plane, with the G2A galaxy being at more than 10cr away. 

The slope of the total mass density profile measured 
by strong lensing provides a value close to isothermal 7 = 
-2.078 (ptot oc r 7 ) at z = 0.2 for ETGs (Koopmans et al. 
2009; Auger et al. 2010; Barnabe et al. 2011). We measure 
the slope in our simulations for the average total density 
within a radius r such that ptot = M t ot(< r)/(4/37rr 3 ). We 
use two points to measure the average slope, at r e ff/2 and 
r e ff. The first striking result is that the measured slopes are 
quite similar to runs with or without AGN, and are close to 
isothermal (see Fig. 17). At z = 0, galaxies without AGN 
feedback show an average slope close to 7 = —2.29 =b 0.07, 
while for the galaxies with AGN feedback profiles are flatter 
with 7 = —1.89 =b 0.16. Note that the measured scatter in 
the AGN feedback case is in good agreement with that found 
in the observations. The difference in the value of the slope 
already happens at z = 2 1 Thus, AGN feedback at z = 2 
has been able to significantly change the total distribution 
of matter. Note that the observational points stand between 
the values of the runs with and without AGN feedback, and 
that the error bars are compatible with both values. 

Combining strong lensing data at higher red- 
shift, Ruff et al. (2011) and Bolton et al. (2012) provide an 
estimation of the variation of the slope 7 with redshift. They 
find that the slope is steeper at low redshift than it is at high 
redshift, with a variation of 0.25 in Ruff et al. (2011) and 
0.60 in Bolton et al. (2012) with redshift. Our simulations 
show such a mild evolution in the slope of the total den- 
sity profile with redshift, and our measured averaged slope 
at high redshift z c± 0.5 — 1 is in better agreement with 
observations for the AGN feedback runs. 

Likewise, when comparing in Fig. 18, the projected DM 
fraction within half the effective radius at z ~ 0.2 in the 

1 Above this redshift the value of r e ff/2 is too close to the reso- 
lution limit to be conclusive for the mass density slope. 
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SLACS data (Gavazzi et al. 2007; Auger et al. 2010) and in 
our simulations, we find a better agreement for the simula- 
tion with AGN feedback turned on (/dm ~ 0.5) than without 
(/dm ~ 0.2), mostly because the stellar mass (or equivalently 
the velocity dispersion) without AGN feedback is too large. 
In the former case, the dark matter fraction is slightly more 
consistent with observations assuming a Chabrier IMF. 



4 COMPARISON WITH PREVIOUS WORK 

A key point of our work is that, in the absence of AGN 
feedback, the galaxies always exhibit a late- type morphol- 
ogy with a large disc of stars as can be observed in 
the synthetic images in Fig. 4. We find similar results 
to Lackner et al. (2012) who use the enzo (AMR code) to 
simulate some overdense and underdense regions of the uni- 
verse: their stellar conversion efficiency is high / CO nv =0.6, 
and their most massive galaxies M s > 10 12 M© are dom- 
inated by in situ star formation /i ns itu = 60 per cent at 
z = 0. They use similar physics (gas cooling, low efficiency 
star formation, supernovae feedback and metal cooling) as 
for our simulations without AGN feedback. In the work 
from Oser et al. (2010) (see also Naab et al. 2007, 2009; 
Oser et al. 2012; Johansson et al. 2012), they employ the 
gadget code (Springel 2005) based on SPH, to study the 
formation of elliptical galaxies in group environments. The 
same basic physical ingredients as for our simulations with- 
out AGN, and in Lackner et al. (2012) are used in their 
work, except they use a larger value of the star formation effi- 
ciency e* ~ 0.1. However the results of Oser et al. (2010) are 
radically different from ours (and thus from Lackner et al. 
2012): they naturally obtain elliptical galaxies with low stel- 
lar conversion efficiency 0.1 < / CO nv < 0.45 respectively for 
high mass and low mass halos in the mass range 10 12 < 
Mh < 2.5 x 10 13 M©, and the growth of galaxies is largely 
dominated by the accretion of satellites with /insitu = 20 
per cent at z = for the halo mass range corresponding to 
ours. Note that, in SPH simulations of M h = 3 x 10 12 M© 
halo by Khalatyan et al. (2008) using SPH, they find little 
difference in their galaxy mass and morphology between the 
AGN case and the no- AGN case at z = 0. The galaxy is 
more red with AGN feedback, and has a larger value of the 
Sersic index, but the galaxy is an elliptical in both cases and 
has the same low stellar conversion efficiency / CO nv — 0.2 at 
z = 0. 

There are several potential interpretations that can ex- 
plain this difference in the results. The first one is that the 
high star formation efficiency adopted can account for the 
difference seen in the accreted versus in situ star formation. 
Higher efficiencies tend to form stars at higher redshift by 
rapidly converting gas into stars, thus, the fraction of ac- 
creted mass should increase at low redshift. We make that 
simple test by running the Gl and G2 simulations with a 
e* = 0.1 efficiency (still without AGN feedback), and we 
found that the accreted fraction of the stellar mass has 
increased. In line with the analysis performed by the de- 
composition of in situ versus accreted origin of stars, we 
indeed find that the galaxy becomes more spherical (sup- 
ported by the velocity dispersion of stars as opposed to 
rotation). However, the stellar mass obtained at z = is 
the same (M s = 5.4 x 10 n M©, / con v = 0.76 for Gl; and 



M s = 7.0 x 10 11 M©, / C onv = 0.66 for G2) as for the low star 
formation efficiency case (M s = 5.6 x 10 11 M©, / C onv = 0.79 
for Gl; and M s = 7.3 x 10 11 M©, / CO nv = 0.70 for G2), and 
cannot account for the difference observed with Oser et al. 
(2010) (/conv ^ 0.2 for halo masses M h ~6x 10 12 M©). We 
observe that both DM mass resolution, and spatial resolu- 
tion are of comparable values in Oser et al. (2010) and here, 
even though it is difficult to compare a smoothing length 
with a minimum cell size. 

Another much more worrying aspect is the different 
behavior of AMR and SPH codes as already underlined 
by Scannapieco et al. (2012) for Milky- Way halos, and then 
by Vogelsberger et al. (2012) who compared results from the 
AREPO code (moving mesh Springel 2010) with gadget in a 
statistical sense. With relatively similar physics, AMR codes 
(or moving mesh codes) tend to naturally produce overly 
massive disc-like galaxies, while SPH codes predict less mas- 
sive galaxies that have already early-type shapes. This is 
in line with our findings: the "usual" modelling of galaxy 
formation (i.e. not accounting for AGN feedback) with the 
RAMSES code leads to overly massive disc-like galaxies even 
in the centres of massive groups of galaxies. This is poten- 
tially driven by a more accurate treatment of the gas insta- 
bilities in grid-based codes as demonstrated by Agertz et al. 
(2007). It can have critical consequences on the amount of 
gas mixing between satellites (cold phase) and the hot ambi- 
ent medium in massive halos which can lead to more efficient 
gas cooling (Sijacki et al. 2012). Moreover, the treatment of 
blast wave explosions (SN or AGN) with thin discontinu- 
ities can have problematic behaviors if the adaptive time- 
stepping becomes too aggressive (Durier & Dalla Vecchia 
2012). 

A careful inspection of the numerical treatment of the 
hydrodynamics and its consequences for the formation of 
massive galaxies remains to be conducted. This is especially 
important as it can lead to an underestimate of the role 
played by AGN feedback in shaping the masses, colours, and 
morphologies of massive ETGs. We reserve this numerical 
comparison for future work. 



5 CONCLUSION AND DISCUSSION 

By means of high resolution hydrodynamical cosmological 
simulations, we have studied the impact of AGN feedback on 
the formation and evolution of massive galaxies in the centre 
of groups of galaxies. The main finding of this work is that 
we have demonstrated that AGN feedback can turn blue 
massive late- type disc galaxies into red early- type galaxies. 
This process happens because AGN feedback reduces the in 
situ star formation in the central most massive galaxy, and 
turns it into a stellar accretion-dominated system. In more 
detail, we find that: 

• The conversion efficiency of baryons into stars in the 
central galaxy is reduced by a factor 7 from / CO nv = 0.7 to 
/conv = 0.1 at z = due to AGN feedback. 

• Galaxies are intrinsically redder with AGN feedback 
than without, with, respectively, g — r — 0.65 and g — r — 
0.45. 

• The DM mass fraction is increased (from (/dm) = 0.53 
to (/dm) = 0.83) with AGN at z = 0, and stellar (from 
(/ s ) = 0.42 to (/ s ) = 0.15) and gas (from (/ g ) = 0.04 to 
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(fg) — 0.016) fractions are reduced within the core of halos 
(in O.lrvir) at z = by the presence of AGN feedback. 

• Central galaxies are rotation-dominated (discs with 
v/a > 1) without AGN feedback, and dispersion dominated 
(ellipticals with v/a < 1) with AGN feedback. 

• The fraction of stellar mass in the halo locked into satel- 
lites is larger if AGN feedback is present. It leads to the 
formation of central ETGs dominated by accretion of stars 
rather than by in situ star formation at z = 0. 

• Galaxies are dominated by in situ star formation at high 
redshift because they are gas-rich systems that rapidly con- 
vert gas into stars. As their reservoir of gas decreases with 
time and the star formation rate is quenched, the proportion 
of accreted mass increases. 

• Accreted stars are deposited far away from the central 
galaxy, while in situ formed stars are close to the centre 
of the halo. As a consequence, the effective radius of the 
accreted distribution of stars is larger than that of the in 
situ component. 

• Because of the increased fraction of accreted material in 
the AGN simulations, the effective radius of central galaxies 
is larger, and galaxies are less compact. 

• The presence of AGN feedback provides a better fit 
to the observed dynamical scaling relations (size-mass, 
velocity-mass, M s -plane) than without AGN feedback. 

• With AGN feedback, galaxies are smaller at high red- 
shift at constant mass, while they are larger without AGN 
feedback. 

• The total mass density profiles in galaxies are affected 
by AGN feedback, and they get shallower. 

We note finally that the conclusions of this work could 
be mitigated by some possible artefacts of our methodology. 
We employ a finite cell size resolution of 0.5 kpc, which can 
possibly affect the compactness of the simulated galaxies 
in particular at high redshift, where their sizes approaches 
the resolution limit. The stellar-halo mass relation at z — 
is compatible with a Salpeter IMF for the most massive 
galaxies and with a Kroupa (or Chabrier) IMF for the least 
massive galaxies. In contrast, the size-mass relation favours 
a Kroupa IMF for the most massive galaxies, while the 
two least massive ones are excluded from the observational 
constraints (Kroupa or Salpeter). This leads to an incon- 
sistency that can be the outcome of a lack of resolution, 
or wrong/missing modelling of the galactic physics, that 
remains to be explored further. As the multiphase struc- 
ture of the gas is not resolved, the formation of massive 
clumps of gas produced by the disc instabilities at high red- 
shift is not triggered (Agertz et al. 2009; Dekel et al. 2009) 
and any morphological quenching associated with this phe- 
nomenon is missed (Martig et al. 2009). However, the in- 
crease in resolution is not enough to prevent the galaxies 
from being too massive and strong feedback processes are 
still required (Dubois et al. 2013). Also, because of finite 
computational time, we limited our sample to six objects 
over one dynamical range of halo mass. More halos span- 
ning a broader mass range have to be simulated to obtain 
better statistical significance on our preliminary theoretical 
predictions both in terms of mean average trends and intrin- 
sic scatter. We stress that, without AGN feedback, none of 
the observational constraints can be reproduced regardless 
of the choice of the IMF. 



We intentionally did not discuss the singular role of mi- 
nor versus major mergers, or wet versus dry mergers, but in- 
stead we emphasised the particular effects of AGN feedback 
in transforming in-situ dominated galaxies into accretion- 
dominated systems. As already highlighted by many au- 
thors, the role of the mergers of central galaxies with their 
satellites is crucial to understanding the size and velocity 
dispersion of galaxies and their evolution with time. It is 
beyond the scope of this paper to investigate the nature of 
galaxy mergers and possible correlations with feedback, but 
we in future work we will perform a dedicated study of these 
effects by employing a more statistically significant sample 
of galaxies. 

Finally, when comparing our results, without AGN feed- 
back, to the existing literature employing hydrodynamical 
cosmological simulations and the formation of ETGs (mostly 
without using AGN feedback), we found our results to be 
consistent with galaxies simulated with an AMR technique 
(blue massive discs dominated by in situ star formation), 
and inconsistent with those obtained with SPH (red ellipti- 
cals dominated by accretion). This conflict calls for a thor- 
ough analysis of the respective merits of the different nu- 
merical techniques and their consequences for the formation 
and evolution of massive galaxies. 
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